Geometry of intensive scalar dissipation events in turbulence 
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Maxima of the scalar dissipation rate in turbulence appear in form of sheets and correspond to the 
potentially most intensive scalar mixing events. Their cross-section extension determines a locally 
varying diffusion scale of the mixing process and extends the classical Batchelor picture of one mean 
diffusion scale. The distribution of the local diffusion scales is analysed for different Reynolds and 
Schmidt numbers with a fast multiscale technique applied to very high-resolution simulation data. 
The scales take always values across the whole Batchelor range and beyond. Furthermore, their 
distribution is traced back to the distribution of the contractive short-time Lyapunov exponent of 
the flow. 

PACS numbers: 47.54.-r, 02.70.-c, 07.05. Pj 



When a scalar concentration field 6'(x, t) is transported 
in a turbulent flow very large scalar gradients are gener- 
ated which can be associated with potentially intensive 
mixing Jj] . Frequently, such large-amplitude gradient re- 
gions exist across scales that are finer than the smallest 
turbulent eddies, a case which is known as the Batchelor 
regime of scalar mixing ^] . Their cross-section is usually 
estimated by a single mean diffusion scale, the Batchelor 
scale 77b, that equilibrates advection by flow and scalar 
diffusion. However, scalar gradients are known to fluc- 
tuate strongly in turbulent mixing. These fluctuations 
are caused by the fluctuating scalar amplitudes and by 
the varying spatial sections across which the scalar dif- 
ferences are built up. Both aspects will cause the strong 
spatial variability of potentially intensive mixing regions. 
Thus, a whole range of local diffusion scales Id, which 
quantifies exactly this variability, will exist around the 
mean diffusion scale t/b ■ Their distribution is interesting 
for several reasons. It can enter mixing efficiency mea- 
sures . The scales prescribe the extension of chemically 
reactive layers in which the combustion of fuel takes place 
Q or the variations of the fluorescence signal as used for 
the measurement of zooplankton patchiness in the ocean 
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In this Letter, we want to calculate this local diffu- 
sion scale distribution p{ld)- It arises from the competi- 
tion of two dynamic processes. On one hand, the scale 
distribution will be affected by molecular diffusion that 
causes a diminishing of existing steep gradients as well as 
the completion of their formation by reconnection Q ■ 
On the other hand, the scales will be determined by the 
statistics of the local advection flow patterns that pile 
up scalar differences. While the strongest scalar gradi- 
ents were related to instantaneous velocity gradients in 
Ref. Q , we will go one step further in the second part of 
the Letter and relate the local diffusion scale distribution 
to the distribution of Lagrangian contraction rates of the 
flow, as given by the smallest of the three finite-time Lya- 
punov exponents along Lagrangian trajectories. The La- 
grangian approach incorporates the temporal evolution 



of the persistent flow patterns that eventually generate 
the strongest scalar gradients or the finest local dissipa- 
tion scales in a finite time. The analysis is conducted for 
passive scalars in Navier-Stokes turbulence, a case that 
is not accessible to analytical treatment. Our investi- 
gations are therefore based on direct numerical simula- 
tions (DNS). In order to discuss trends with Reynolds 
and Schmidt number, we analyze five different data sets. 
Very high-resolution calculations (with up to 1024^ grid 
points) have to be conducted with resolution constraints 
that exceed the usually applied ones such that the finest 
scales of the turbulent mixing process are resolved prop- 
erly. The fluctuations of the flow are sustained statisti- 
cally stationary by a large-scale random forcing and the 
passive scalar is driven by a constant mean scalar gra- 
dient. We consider the scalar dissipation field Q which 
probes the magnitude of scalar gradients and which is 
defined as 



ee(x,i) = K|V6l(x,i)|2 



(1) 



Here, n is the scalar diffusivity. Experimental studies on 
the geometry of scalar dissipation fields are very chal- 
lenging since gradients have to be measured and only a 
few exist 

Figure 1 shows a two-dimensional (2D) slice cut 
through a DNS snapshot of eo (x, t) . We observe strongly 
folded filaments. The points that form the level set of 
largest dissipation amplitudes 



{x : ee > C{ee) or |V0| > ^C{\Ve\^)} , (2) 



are redrawn in red. C is a real constant. The resulting 
filaments are cross-sections of thin sheets in which the 
maxima of scalar dissipation are arranged in the three- 
dimensional volume Q. A closer inspection of Fig. 
unravels various length and thickness scales of the fila- 
ments. The filaments arc curved and tightly clustered 
in certain locations thus posing a challenge of separating 
each curved filament and computing its accurate varia- 
tion scales. 
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FIG. 1: (color online) Contour plot of a two-dimensional slice 
cut through the instantaneous three-dimensional scalar dis- 
sipation rate field. Level set Lc for C = 4 (see Eq. 1^}) 
is replotted in red. Data are from a pseudospectral simula- 
tion of the advect ion- diffusion equation for the passive scalar 
in combination with Navier-Stokes equations for a statisti- 
cally stationary, homogeneous isotropic flow at a resolution 
of 1024^ grid points in a periodic box V = (27r)^ Q. The 
Schmidt number is Sc = v/k = 32 and the Taylor microscale 
Reynolds number is R\ = -y/ 15/ {^{e)) {u^) = 24 with (e) 
being the mean energy dissipation rate and v the kinematic 
viscosity. The Batchelor scale r]B = •q/^fSc is resolved with 
2 grid cells. The viscous Kolmogorov length rj = i'^''^/(e)^''* 
measures then 11 grid cells. This spectral resolution is larger 
by a factor of 4 than the one usually adopted. 



This is done here by a fast multiscale clustering al- 
gorithm 12J, which is based on the Segmentation by 
Weighted Aggregation (SWA) algorithm motivated 
by Algebraic Multigrid (AMG) [ll]. The algorithm as- 
signs data points into clusters starting at the finest res- 
olution level, the grid spacing. All points e Lc are 
therefore gathered in a so-called proximity graph that 
contains their location and their connectivity to other 
points of Lc as quantified by the inter-point weights 
which probe the nearest graph neighbors of each point 
only. Since the information about proximity is kept as 
one moves from finer to coarser resolution level, one 
can perform a fast recursive principal component anal- 
ysis (PC A). The resulting eigenvalues characterize the 
length and width of the point clusters. Additionally, 
strongly curved filaments have to be decomposed into 
sub-filaments by applying a local convexity criterion 
along the filament. Figure [3 shows the reconstruction 
of the filaments from Fig. 1 and their division into sub- 
filaments. The local dissipation filament thickness. Id, 
is then given by nothing else but the smaller eigenvalue 



FIG. 2: (color online) Reconstruction of the red colored fil- 
aments as shown in Fig. 1 by means of the fast multiscale 
clustering algorithm. Long filaments are composed of several 
subfilaments that are colored differently. 




FIG. 3: (color online) Distribution of the local cross-section 
thickness Id of the scalar dissipation rate filaments for eg > 
4(ee). Left panel: Probability density function (PDF) 
p(ld/riB) for three different Schmidt numbers at R\ — 24. 
The dashed line corresponds with the theoretical value of the 
most probable scale 1^ — ywi/jAal which will be discussed 
later in the text. Right panel: PDF pild/ff) for two different 
Reynolds numbers at Sc= 32. 



which follows from the PC A of each sub-filament. 

FigureOlshows the probability density functions (PDF) 
of the local filament thickness Id for different Schmidt 
{Sc) and Reynolds numbers {Rx). The distribution will 
depend on the cut-off level C, but the physical picture 
will not change since C is fixed with respect to the mean 
scalar dissipation rate throughout the analysis. Local 
thickness values within the whole Batchelor range be- 
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FIG. 4: (color online) Probability density function of A3(t) 
for different times. The data are for Rx — 24. The variable 
z^{Xs{t)^{Ht)y)/as{t) withaaW = ^iiMW^I^Mf) 
is chosen for comparison with a Gaussian distribution. We 
found A3 ~ -0.425 and cr|(t) = a^^/t with ~ 0.44 for t > 
jAal^^. The distribution is evaluated by 2.5 x 10^ Lagrangian 
tracers that are initially seeded uniformly. The distribution 
for R\ = 10 behaves qualitatively the same. 



tween ijb and 77 and beyond are found, indicating that 
dissipation maxima are also related to scalar gradients 
across inertial range scales. The 2D analysis does not ac- 
count for the spatial orientation of the sheets with respect 
to the cutting plane. As demonstrated in this will 
affect only the tail for large Id- The left panel of Fig. O 
compares the PDFs for three different Schmidt numbers 
at a fixed Reynolds number R\ — 24. With increasing 
Schmidt number stronger jumps of the scalar concentra- 
tion across finer thickness scales become more probable 
since diffusion is less dominant. Consequently, the most 
probable thickness Imax-, i-6- the maximum of the distri- 
bution, is shifted to smaller values, but remains always 
larger than the corresponding Batchelor scale. This sug- 
gests that the formation of so-called mature scalar gra- 
dient fronts with a thickness ~ 773 is a subdominant pro- 
cess. We see that the three PDFs overlap when rescaled 
with TjB which implies that Imax ~ Sc~^^'^ on the basis of 
our data and even though Sc = 2 has no real Batchelor 
range. In the right panel of Fig. O we compare distribu- 
tions for two different Reynolds numbers at fixed Schmidt 
number. Both PDFs rescaled with the corresponding rj 
overlap again to a large fraction, except for the very fine 
scales. Their higher probability with increasing Rx in- 
dicates a more efficient stirring at the smallest scales. 
Our data show that the scale Imax follows now the same 
dependence with Reynolds number as the Kolmogorov 
scale, i.e. Imax ~ R\ ^^'^ This result is similar to finding 
in fo'^ ~ 1 and does not change for Sc^ 1. 

The filament thickness distribution is now related 



to the distribution of finite-time Lyapunov exponents 
(FTLE) Xi{t). They measure the local separation be- 
tween two initially infinitesimally close fluid elements 
along the Lagrangian trajectory. Therefore an orthog- 
onal frame is attached to each tracer and the three sep- 
aration vectors evolve as d(5r j*' (t) /dt — Ujk (t) Sr^j^^ (t) 
for j, /c = x,y,z and i = 1,2,3. o'jkit) is the rate 
of strain tensor along the Lagrangian trajectories. The 
FTLEs follow to Xi{t)=l/t log(|5rW(t)|/|<5rW(0)|) where 
the Gram-Schmid method is applied consecutively to 
three Sr^-^^t) (see Appendix C.3 of The exponents 

will vary from trajectory to trajectory resulting in sta- 
tistical distributions. Their means are the global FTLE, 
(Ai(t)), which converge for long times to three numbers, 
the asymptotic Lyapunov exponents A^ = limt^oo (Ai (t)) . 
Due to incompressibility, J^^^i^ii^) — 0. Our inter- 
est is in the formation of thin dissipation (or gradient) 
sheets where expansion in two directions is present, i.e. 
Ai(i) > and X2{t) > 0, and contraction in the third 
one, A3(t) < 0. The distribution of contraction rates 
that pile up scalar gradient maxima and form the distri- 
bution of the local diffusion scales follows consequently 
from the PDF p{X3{t)) as shown in Fig.0] We see that 
the cores of the distributions for different times collapse 
to a Gaussian profile within ±2a3(t). Since the standard 
deviation 0-3 (i) — (Toc/Vi for larger times the contraction 
rate X^lt) will get more and more concentrated about 
A3 as the time advances. Base d on A3 the most proba- 
ble thickness is given as l*^ ~ ^k/|A3| [l^ which arises 
by equilibrating contractive strain and diffusion. For all 
data analysed this scale is at about the maximum of the 
thickness distribution, ~ I mn.x (see left panel of Fig. 3). 

Hu and Pierrehumbert [l3| pointed out that the 
asymptotic A3 alone is not sufficient to explain the for- 
mation of the finest diffusion scales when the flow is time- 
correlated as being the case for the present studies. The 
time that a scalar blob or filament experiences a persis- 
tent strain pattern along the trajectory can be estimated 
by a characteristic decorrelation time which is given by 
Tc ^ ^3^- For larger periods the persistent compression 
of the blob will disappear or the sheets undergo diffusive 
destruction or merging. This local scenario will appear 
repeatedly and causes a stationary thickness scale distri- 
bution. Consequently, we have to take into account the 
evolution of distribution of the FTLEs over such periods 
in order to study the formation of the gradient (or dissi- 
pation) sheets. The distribution of the sheets is related 
to the characteristic distribution of the FTLE over times 
i < by 
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dlop{lo)g{\og{ld/lo)/Tc), (3) 



where the distribution of the FTLE is approximated 
by a Gaussian g{X3,t) = (1/^2^3(^)2) exp[-(A3(i) - 
(A3(i)))^/(2cr3(t)2)] andp(Zo) the (unknown) distribution 
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FIG. 5: (color online) Distribution of the local cross-section 
thickness Id as reconstructed from the distribution of the 
FTLE (?(A3) via Eq. for different times. Data are for 
R\ = 24 and Sc = 32 (see Fig. |3J . Two initial distributions 
are shown: (a) uniform distribution p(io) = l/{2r]K — Vb)', 
(b) delta function p{lo) = 5{lo - 1.81*^) (see Fig. 3). The grid 
spacing A and Kolmogorov scale rj are indicated as vertical 
arrows. 



of the initial thickness scales. The initial distributions 



p{lo) taken here mark both ends of a spectrum of possi- 
ble choices, the localized case given by a delta- function 
around the maximum of p{ld) and a uniform distribution 
of initial scales over an interval, respectively. Both cases 
can be carried out analytically and result in scale distri- 
butions of log-normal type. Figure 5 shows the result- 
ing PDFs for Rx = 24 and Sc = 32. The distributions 
agree well with the data for most Id and for time periods 
t ^ 1^3 1 The uniform case fits slightly better for the 
left tail since sufficiently small scales are present initially 
that can be steepened further to smaller cross sections. 
The strain which is accumulated over short times only 
seems to explain the formation of the finest scalar gradi- 
ent sheets. Nearly, the same time scales and distributions 
were found for the other Schmidt and Reynolds number 
values. 

In conclusion, we have determined the distribution 
of the cross-section extensions of the scalar dissipation 
maxima in Navier-Stokes turbulence. These scales corre- 
spond with local diffusion scales and take values across 
the whole Batchelor range and beyond. By means of the 
distribution of smallest short-time Lyapunov exponent, 
the diffusion scale distribution can be reconstructed. 
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